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ABSTRACT 


The military officer personnel system is modeled as a 
semi-Markov process. The probability distributions de- 
Seri vimte=starce transitions are=developed in terms of policy 
Vetrlavees. ““H@rrects Of chanees ain promotion policy and force 
size are studied for both a closed, steady state system and 
an open, transient situation. Moments of the joint distri- 
bution of officers in each rank are presented for all cases. 
A method for approximation of transient results from known 


iste rmmlc1ons 1s formulated. 





TABLE OF CONTENTS 


iene Wee Net = ee = 4 
Mi MOMEINEORMUWATION --=9=4e5=------ 22522-22222 s222- 6 
III. STEADY STATE RESULTS ----------------------------- 15 
IV. TRANSIENT RESULTS -------------------------------- 27 
LIST OF REFERENCES --------------------- wenn eo -------- 36 
INITIAL DISTRIBUTION LIST ----------------------------- 37 
FORM DD 1473 ---------------------- eo ee eee eee 38 





I. INTRODUCTION 


Personnel systems composed of a formal rank or grade 
Senueture that utilize a prenetion procedure to move™per - 
sons through or out of the system are commonplace in modern 
eerae Orpena zations. In particular, the Officer Corps of 
cweimomeche=mulaitary services*provides a classic example 
Ot Bhis type. 

lmineary officer promotions are controlied by federal 
lawsssuch as the Officer Personnel Act of 1947 and the 
Officer Grade Limitation Act of 1954. Additional constraints 
may be imposed by service Secretaries, frequently acting as 
approving authority for the respective service's internally 
@eniem@ated promotion pian. 

Seteleormye ceilings, ~otten fureher reduced by service 
Polvetes, CX1St on the maximum number of officers holding a 
particular rank. Promotion boards meet annually to select 
Siieers £Or promotion to the next higher rank, Total num- 
Bers to be promoted are established in accordance with pre- 
[mosmend projected attrition. Ihe size of the prometion 
meme —semese officers to be considered for selection ~- is 
based on a promotion opportunity Dareneter . 7c 20 es 0.0 Dr On 
motions might be obtained by a 50% opportunity from a 1000 
ime zonemeor by a 80% opportunity from a zone of 625. 

Time in grade - the number of years an officer has 
Spent 9 1n mis current rank - is a pecondany consideration in 


fieme the size of the zone. That is, the may#mum speed 





with which an individual may move through the system is a 
function of the size of the zone each year.!? 

Pus, there exist three policy variables, (1) total 
number in each rank, (2) promotion opportunity to next 
rank, and (3) time in grade, which control the functioning 
of the promotion system. 

The purpose of this paper is to present a mathematical 
medel of the flow of personnel through this type multi- 
grade system and to study both long range effects of a given 
promotion policy and short term effects of a change in 
policy. The paper is written in four parts. In Part II 
we formulate a semi-Markov process model and define its 


variables and properties. Part III contains steady state 
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demuve a method to approximate them. 


1 
Accelerated promotions "below the zoneé'' do occur. They 


account for less than 5% of the total promoted but will be 
Miciwaecdmin all probability distributions disc isséd in this 
paper . 





II. MODEL FORMULATION 


Personnel flow models in which movement from state-to- 
State in successive time periods is assumed to follow a 
Markov chain have received considerable attention in the 
literature (e.g., Bartholomew [1], Gani [3]). An extension 
of a Markov chain model which accounts for some of the non- 
Markov properties of personnel systems is to model such a 
System as a semi-Markov process (S.M.P.), 1.e€., a process 
whose sequence of states forms an embedded Markov chain at 
the transition points, where the system changes state, and 
in which the lengths of time between these transitions are 
random variables (see Pyke [6]). 


a] a i 
s 46 


Eet us fixss@® tormulate an “open model,” Ynat rast 


> 35 ay te 


pele 


ace sj.) 3 = 1,2,..-,6; 1s the rank of an active military 
officer and state 7 is an absorbing state containing those 
Giercersewho are retired, resigned, deceased or otherwise 
[Orerore tie service ,"' then ranks 1 through 6 are™a set or 
transient states that form the "system.'' Persons leaving 
this set of states must be replaced by new recruits into 
ile. Swen aut asenotetomempty over tame. Figure 1 il- 
lustrates movement in an open model. 

Byecemtirolling the input flow of new personnel we cam 
control the total number in the system. If that control is 
such that in any time period the number of departures is 
equal to the number of arrivals, the number in the system 


remains constant. Since in this paper we beg 1 by studying 





a "sy oevem Of fixed size fin Part “1V ths assumption is re- 
axed) we can "close" (see Cinlar [2]) the model by arti- 
ficially closing the feedback loop shown as a dotted line 


jr rreure 1. 


NEW 





Figure 1. Personnel Movement in Open Model. 


Now consider a system with 7 states, where state 7 is 
artificial and represents the absorbing state for those 
leaving active service and from which a new recruit comes 
COmemeer tine active States and wéplace the departins officer 
(although he does not necessarily enter the state from which 
Clemiecpameire OCCUrTed) . 

As we shall see, the times taken to move from state to 
state in the model play an important Pole, and for st@ttes i 
firomemeo these times™are dependent on»=polacy decisions = by 
fixing the time to move from the artificial state 7 back to 
the active states at zero (instantaneous replacement), we 
Gaieeeeieovemstate 7 from the system entirely and close the 


HGGpein Fieure 1 without it. Thus We are lef With a 





Glesead System of 6 recurrent, non-null states in our semi- 
Markov model. 

itjees a system of this type which will be considered in 
this paper. We shall allow inputs to the system to occur 
only at rank 1. In Part IV we evaluate increases of overall 
size by new cohort inputs to the fixed size. 

Weenowepursue the system containing» states 1, e2j059e4, 
5, 6, the officer ranks from 01 to 06.* As time passes an 
individual officer moves through the system from state to 
State. A transition occurs when he moves to the next higher 
state via promotion or when he leaves the service, in which 
case he returns to state 1 and is "replaced" as previously 
described. We now assume that the sequence consisting of 
pemeseermociacely rollowing transition” pornts frags a 
Mapwomsehaan. ltuwis significant toonote here thatthe wa- 
lidity of this Markovian assumption is established by the 
Special structure of the mrlitary promotion system * an 
Officer either moves up to the next rank or he leaves the 
system. Thus knowledge of only the present rank provides 
complete information concerning the probabilities for the 
Meets Uramscition. 

Times between transitions are random variables that de- 


pend only on the current state and the state to which the 


*Rank/Grade USA-USAF-USMC TITLE USN TITLE 
O01 Second Lieutenant Ensign 
02 First Lieutenant Lieutenant. jace 
WS Captain Licutenant 
04 Major Li utenant 
| Cummander 
05 Lieutenant Colonel Commander 
06 Colonel Captain 
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transition will be made. Implicit here is the assumption 
that accelerated or delayed promotions do not affect later 
behavior. 

Fomiabiy henge considemmthemscate space {1, 24. 35) 4.05" 
Of, thewtfirst six officer ranks. Let {Z, (t) = i} be the 
Ewent Officer k is in rank i after a time t in service. 
Then the stochastic process {2, (t), t > 0} is a semi-Markov 
process (S.M.P.). The stochastic process that records the 
number of times the S.M.P. has entered each of the system 
States) 15 Calléd a Markov Renewal Process (M.R.P.),) as nomed 
Dywueees [6.7]. The S.M.P. and the associated M.R.P. areguene 
ally discussed together in the literature. However, the 
Special structure of a personnel system is such that the 
POubLOoWwine GiSCUSSICN Willi be iimited to tme S.M.P. 

We shall assume that each officer in the system fol- 
lows an independent path from the same S.M.P.* If we let 
X5 (t) be the number of officers in rank j at a time t, then 
{X(t), t > 0} where X(t) = (X,(t), X,(t),...,X(t)) , is 
Gamiule¢il-particle semi-Markov process penerated by he =. = 


We can now define 


process next makes a transi- process has 
Q. . (t) = Pp LON MEO | @dndecitsTOCcuUrs Just Cntene 
J iia tine << ¢ i a 
(1) 
Gttiecer k N€xt oes sto ran. he has juse 
eee lj wetar this OCCURS Siimagetame attained rank 
St 1 


7The state space could, of course, be des-ribed in gen- 
eial. it 1S limited to the specific case to | discussed £0 
Simplify subsequent notation. 

*See Heyman [4] definition of a Multi-particle SMP. 
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H, (t) 


I 


State at Citwen t alien] 


officer's *rankweat time 
Cis j 


ce) = 3 (Zao) =a 


areransition Wililieeoccur 
Die alee amouUn C oOnemesl Cama 


Gukiecer' Ss Lancy wael 
change in an amount of 
enc < t 


next rans Perens oceirs 
in an anoume O14 time 
< t : 


Cime, Obllecerespecias=im 
rank i<t 


first transition Inte 
State } Occurs In a Cime 
<a 


Drecess 
Started in 
Stace 1 at 
time 0 


officer 
started in 
Task ~ dee 
time 0 


| 


process has 
just entered 


i and will 
NeXt eiaee nce 


(3) 
off1cer Mas 
just attawmed 
rank i anid 
will next 
noid em = 


kL. 


process has 
VUSt Eller 
i 


4) 


( 
he has just | 
attained { 


rank i 


process 
Starts -limue 


(5) 


Tt sshoutasbe noted that 655) does not have a meaningful 


explanation in terms of the functioning of the officer sys- 


Com. Wowever , 


Fesm@utcmiicuitseits inclusion ingthe development. 


1tseusetulness se thewecri Va balon 
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Samneec Qi 4 igmthe prepabusity that the next trams. - 
tion is into state j given the process has just entered i, 


we define: 


Qi 5 (2) Ge Pij 
so that PG is the fraction of those who reach rank i who 
ever attain rank j. 


Denoting the matrix of Q; 5 (8) by @(t) [we will tence - 


forth use the notation ay = (En, ve write 


© 144, Pao 2 0 0 0 
ee 060 Pog . : 
S A324 0 0 Pra 0 0 
. 7 ee er 
dei) aes 0 Y 0 P56 | 
= 14 0 0 0 0 0 


where q34 7 i} is the fraction of those who reach 


~ Piitl 
rank i but are never promoted to rank itl and are replaced 
by new officers in rank 1. 

It should be re-emphasized that this matrix defines a 
closed system - a system in which the number of Officers. 
remains fixed with no extra entering or leaving the system 
over time. Thus, mathematically we have a certain tract Lom 
leaving each state and returning to state Tewath the imeenw. 
pretation of instantaneous replacement by a new officer. 

Recalling that oF ee is the joint proba-ility of the 


next teetisition from i to j and in a time les» than t We have 
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Sainrce ene is the probability that the next transi- 
tion is into state j given the process has just emuemedm: , 


we define: 

Q3 5 Ga) a2 P3j 
Som ena t Pi; is the fraction of those who reach rank i who 
ever attain rank j. 


Denoting the matrix of Freee, by Q(t) [we will hence- 


forth use the notation Q(t) ~ (Q;5,0t))] we write 


/ 
i ae : : : 
< 491 0 P53 0 0 0 
= dz 0 0 P24 0 0 
a(e) = ao 
| Gay 0 0 0 Pas 0 | 
| 
Gs, °° 0 0 0 P56 
1 0 0 0 0 0 
where q;1 = UES Ret is the fraction of those who reach 


rank i but are never promoted to rank itl and are replaced 
by new officers in ram. 1. 

It should be re-emphasized that this matrix defines a 
closed system - a system in which the number of officers N 
roVieuine fixed myth no extra entering oF leaving the system 
over time. Thus, mathematically we have a certain fraction 
leaving each state and returning to state 1 with the inte we 
pretation of instantaneous replacement by a new Ofticecr. 

Recalling that Oe is the joint proba’ ility of tite 


next transition from i to j and in a time les» than t we have 


————— 


ll 








..(t .(t 
ss all ey : 0 - 
and thus 
6 6 
H(t) = ) py; FR;(t) = D. Q5(t). 
jel j=1 


It can be shown (Pyke [6]) that 


6 le 


Pig lt a OG: ya H, (t) Pex \ i Pig FE XD Q5 4 OD) 
KkK=1 550 


where mae is the Kroneker delta. 


If we let H. (t) = j]- H. (t) and gg We x Bee india. 


cate the convolution we obtain 


6 
Pools) Sees eKCs a Se, Se Se aes Il ee 
J =) 


The matrix form of these equations 1s 
EtG) (= 2 (eye ae eee (7) 
where fel) Loa diaconalemac lx On H, (t), and “Of t) % 9a) 
iomene matrix convolution. ” 


A similar relationship between eg CS and gle) can 


be defined as 


6 ie t 


8 C8 - x J Ca ee -f Hee Bee EN Bars: 
O O 


with the matrix form 


°A convolution of matrix valued functions is similar to 
ieeriomium~wetplication except the product of G€leiientS 1S fee 
placed by the convolution of elements. 


EZ 





G(t) = @ # G(t) + Q(t) - @ * @(t). (9) 


in terms of the policy variables of Part I we cam now 
define (Ross [8]) the average amount of time spent in rank 
i as 


Lee) 


ee a | td H, (t) - f H. (t), 
O 


O 


and the average time spent in rank i given the next transi- 


tion is to rank j 


1% co 
dF “i 
O O J 
romaliy, the mean first passage time from i to j is 
foe) 
f 4- 2 ce f 4. ~ = ne 1 
O 


merowan the concept of first passage is of little direct 


use in this model, it should be noted that G -) is the 


it 
distribution of lifetime in the system. For the most part, 
however, G(t) will be used to derive other results which 
a@yeein terms Of the policy variables. 

In the officer promotion system we can consider 9(t} 
(and thus H(t), Ws and Nig SS known since the Q; 5 0) are 
determined by promotion policies and officer behavior pat- 
terns. Parts III and IV will discuss the unknown P(t), 
whose first row is of primary importance. RS) 


1,2,...,6, gives the expected number of the N officers who 


amen in wank j} at time t. For large t, the Pi. t) cam be 


i 





Wnmempreted as the fraction of officers expectedmto be am 
tamk 3. %Ihus, to evaluate the lonp term effect of a givem 
poprey, the lamilting.féerm of P(tjeas requimredeand itsedeter- 
iMaavatiion directly from (7) is formidabhe. The imtroduction 
of thesdastribution G(-) and in particular (9) wibl simplify 
the task. 

Meascumementswof the lamutangetoemm of P(t) ‘comlid eet 
course, be done directly via a cohort type model (e.g., 
Marshall, [5]), but these direct measurements do not give 


any information as to how changes in policy affect the 


eee 
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Tift. STEADY STATE RESULTS 


Ket Us fiest examine mee output of this model after a 
deme perioa of time; that 1s, in the steady state. 

Ins tiecwsteady state, themilimating formof Pi@e) come aum 
the distribution of the N officers over the 6 states. We 
emall see that its solution is indirectly in terms of pol- 
hey weeetables and will illustrate how the effects of pol- 
ey Changes can be measured in terms of the parameters of 
mie model. 

Dewwnine the Laplage Stieltjes transform (LJS.a1. jen 


jive function F as 


F(s) -{ ede) s > 0 
O 
we solve simultaneously equations (7) and (9). 
Porstemtaking transtoerms Of (7) we obtain 
P(s) = Hp(s) + Q(s) P(s), (10) 
widen cCan also be written 
A =k _ A iad i | 
2 - ee ae Cs 
Likewise, transforms in (9) yield 
G(s) = Q(s) - @(s) Gp(s) + @(s) G(s) (12) 
which can be simplified to 
“aA -] aA “A -] a 
Eras) ) = 2 pees) (Gea n(s)) = (13) 
Potton! )jesand (L3)swe have the result 
P(s) = Hy(s) + G(s) An(s) (14) 


as 





where Ayn(s) 1s a diagonal matrix with elements 








H, (s) 
ae ee 
te analyze the steady state results we now find lim P(t) 
~ C70 
Se cquavalently the lim P(s). Thus, 
$>0 
lim P(s) = lim #,(s) + lim G(s) lim A,fs) = 0 
D D 
s>0 s>0 s>0 s>0 
+ 1 lim Hits) . ee a 
s>0 = 
Craie? 
: aS) 
0 - 6 
| Gog CS) 
pa | 
where 1 is a 6x6 matrix of 1's. 
Now, using L'Hospital's rule we see that 
H. (s) re (a) 
lim =—=—= = im ——— ee 
! 
S>0 G; 5 (s) s>0 G:5(s) Wag 
so that finally we can say ‘i 7 Le 
mal 22 66 
A A u u H 
eee) Solin Pes) Sf ee AW. als) 
t7© s70 Mil ge cor 
i. EZ See ee 
mail = 22 "66 | 
We will now define (see also Ross [8]) th terms in each 


normore en:s matrix of limiting probabilities as 


1G 





Py = Lim P35 (t) -, ae ee toe falta) 
In 
mirat 35, im the steady state the fraction of officers in rapk 
j equals the ratio of expected time in rank j to expected 
Pemiiemence time from rank j bagk to j. 

However, since it was previously noted that, with the 
exception of Wiq the a have no recognizable meaning with- 
in a personnel system structure it is desirable to develop 
an expression for Pep = es which is in terms of policy 
variables. 

Consider the embedded Markov chain and the long run 
Proportion of tranSitions into each rank j; that is, the 


solution vector m., j) =e, 2 EO 


TT = 1 
1 6 j-l 
1+ 3 Pk k+1 
jz2 k=l 
al 
1) Pk k+l 
> ee Z 
T = Teal x LOS; 6 
= ». eye eal 
j=2 =I, 


itmeanealso be shown’ that 


°These results on the expression of P. in terms of 7. 
emcees tc Ross [8], pp. 104-105. J - 





ee oad », Pik “kj 


kf} 
and combining this with Ps = Bede we have 
Te lhe 
ar ee ne es 


J 


6 

T. u. 
uray 
i=l 


Bwstibstrtuting the values of - derived above we obtamn the 


TI 
ia 


Ime y result 


P 
_ k=l “kpkeel ae 
's ee Jo Ly 2 Ommele a 
Ca ». Mi T\ Pk k+l 
1=2 k=] 
Weer ene Numerator cf P, 1s defined as u,. Hence, eaqwa- 


oe 
ny 
oS 


maton wie) rs an expression for the fraction Of Off 1CeTrs im 
miemsy stem that are in rank j in the steady state, in teams 
Gf mean times in grade and elements of Q@(~) - six Ws and 


five P; Thus each i US eae tunet ron  OfeaG. NOs ey eile Ven 


5 
Parameters, but can be computed without use of the detailed 
Seructire of See. On: H.(t). 

Tumse interesting to nete the effect of parameter 
@Saangescman (17). For any j it 1s easy to verify ‘that when 


2 ismbimerecased(decreased) : 


aot) 
e iMercasesdecreases isu 4 


Pe decreases(increased) 1 > j. 


Tov evatitate the sensitivity of 2 to changeesein. agi 


PicIManeiewe let us first consider the case of a change oFea 
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Seeeriuved amount. The denominator of (17) 
BM ee 7 2s MAP poe ose sa PSP TOP oP sara 
OGe 12e os 5A) ees 


mmen 1s affected most by a change in Wy» Since thier coer 
Stents of Hs 1 = 2,..6 jane less thangunity. Simei, Pe 
becomes less sensitive to this type change in Me ess fist 
Ggeases from 1 to 6. However, if a percentage change in 4 
Us 1s considered then the actual values of the us and Ps; 
Tildtaaewergmine the u, to which Ps isemest sensitive. kar wa 
See Oot Parameter values that approximate current policy in 
all military services, U5 1s thesmest critical paramerer 
mor Percentage change analysis. 


The case where all Ws=u and all age provides fumines 


misbeht into the effect of parameter changes. Here (17) be- 


comes . 
j-1 ja ee 
P; = UP = Pr... Pe ; 
6 
=i i 
u + up an 
i=2 


which shows that in this special case = 1S Inde pende muse 
bee Fugiize 2 gives the plot of “3 versus j for various values 
GL ape. 

Returning to the concept of the MSMP(Multiparticle) we 
can now consider the joint distribution of the number of of- 
fi@ems im €ach rank. Since P. is the probability of reaching 


Bpeatemimamasimole trial, and the movement of “ach otiveer 


eS 





IN RANK 


FRACTION 


Fisotiee 2 . 


DeStribitrrenm: in 


Rank for DP; 
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isan independent path, the joamt distribution is multi- 
nomial; that is, X(#) = wai. ee is distributed 
multinomial with parameters N and P where P = Ce 
a row from P(o). 


Thus, 


6 
P| Ay = Xp Xy = Xqo+- kg = Xe » i elt, 


= 
a 


i 
ns. 
ta 
“4 

Re, 
4 
~ 
bo 
u~ 
o> 
On 


For the steady state effect of a given policy we can 


then say that 


Py 
Po 
D | 
; 5 
Eee) |= CN 
4 
Ms 
ES 
and (19) 
P,(1-P,) -P,P, -P|P. “P,P, -P,P, “P,P, 
-P,P, P,(i-P,) -P>Ps “P,P, -P,P. -P,P, 
-P,P, -PzP,  Pz(1-Ps) -P3P, -P =P. “PP, 
V[X(2) ] = Ni. . - e z -Pp Pp 
Se we ey Oe 4P6 
: : : : : -p_P 
PoP, = PSP a PoP. PoP, Pg (1-Pg) 5s 
“PoP, -PgP -P,P. “P,P, -PpP. Py (1-Pg) 


Hence, Xx. and ner are not independent, but are negatively 


Gomrelared such that 


Zi 


2 
a 





V(1-P 5) (1-P5) 


For large N the distribution of X(~) is approximately 
multivariate normal. Thus the distributions of each X. are 
approximately univariate normal with mean and variance given 
by the appropriate elements of the mean vector and covariance 


Matrix above; e€.¢., 


Ss : Oi 0 aie 
/NPZ-Py as 


This enables us to make probability statements about the 
number of officers in any rank in the steady state. 


Returning to the. simplified case of just two parameters, 


we illustrate the correlation between the numbers of officers 


in two different ranks. Namely, for p = .7 
oe 
6 
p(X, Xo) = : 
1-P, 1-P,) 
| ONL 


Tavis, the effect of a change in the number in rank 3 on the 
number in rank 6 (and vice versa) is small. 

im menot the purpose of this paper to describe de 
tailed methods of finding the distributions or estimating 
the parameters used in this model. It must be noted, how- 
ever, that the us and Pi; are not themselves the policy 


variables of a personnel system as described in Part I, but 


mabnciewane functions of those policy variables Tne Pi; 


oz 





repmcseie tie traction of those that Enter rank 1 that ego 

on to rank j. They say nothing about attrition that removes 

Serrecers trom the system (sends them back to rank 1) priceer 

to their consideration for promotion. For example, consider 

P>z- Let us assume that the expiration of tours of obligated 
meiamce lho) £Or a groupeot OLficems occurs during thestime 

they are serving in rank 2. If 2/3 leave the service at 

EOS or before with terminal grade 2, only 1/3 are later con- 

sidered for promotion to rank 3. If the promotion opportu- 


Mey eto rank 3 (for that 1/3) is .9. then we may approximate 


Poz = (.33}(.9) = .297. 


From the discussion in Part II on the various moments 


t canedewseen that 


fie 


used in the model, 


6 
ae > Des naj 1? = <a Ore (20) 
sail 


Tire re deorincmererexpCCled GCIMe= Spent In Tanke 1 Jo1VCm omc 
ext nemesis to be j. For each i the summation in (20) con- 
ees Onuy two terms, 1.e€., Pig Msq + Pag eg47- The Nesey? 
expected time in grade of those to be promoted to the next 
rank, might be considered an.actual policy variable. The 


nN; , are consequences of decisions on the Ns promotion 


sled Oe 
poOraurrty, and the various other forms of attrition. As> 
suming three years in grade 2 before promotion, 1.5 years 

in grade 1, three years obligated service, and Pio = -50 we 


may approximate 


tee = G7 0)1.S)) eras = eacor 
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Vomeecncluge this Section Pemeus “consider ame con- 
plete example of this steady state analysis using a set 
of parameters developed from the above thoughts and rough 


approximations to current military policies. 


bet 
oe lee eg aE 
U5 = 2 P53 = Ae 
Leet a 
Wy = 5 Pas = 5 
Se eC 
eee 


Using (17) we calculate the fractions expected in each rank 
in the steady state. Table I presents these results, the 
Home Gemmeerrect Ci increasing u, to Six years, ane thie 
Seevial bo7. distribution in the U. S. Marine Corps for “cone 
barison. ’ 


TABLE I. FRACTIONS IN RANK 


Approx. Increased Actual 

Rank Parameters ibis USMC 
il 295 .188 ILS 

2 nae ~2 59 . 319 

3 sis 30 se 

4 led Pigs LSS, 

5 PASS | mys .080 

6 .093 094 2055 


Eromm@lojmand the above results 1t tollows: 


7Source of the Marine Corps data is Combined Lineal 
Bist @@f Officers on Active Duty in the Me@rine Corps. 
merronnany 1971, NAVMC P-1005-PD. 
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Inspection of the covariance matrix shows V(X,) ase Cine 
tarcest and V(X6) as the smallest. However, from Table I 
we note that the numbers of officers in ranks 2 and 6 are 
largest and smallest respectively. To provide a more valid 
eEompiarison of the variation in each of the ranks, we employ 


Mie wCOCTLTIGIents Of Variation: 


YNP.(1-P. | 


cv, = nP 


Using N = 20,000° these coefficients were computed using the 
results from the approximate parameters and are listed in 
ape 10, This table also lists for each Me the Ianeervar 
extending two standard deviations on either side of the mean 


in which we expect .955 of abil realizations to fall. 


Thus, we may conclude that when normalized by the num- 
ber of people in that rank, the variation is greatest in 
rank 6 and decreases as j decreases to 2. This is, of course, 
iImiGuittvely appealing since the officers in state 6 have jem 


Pemtenceastne largest number of state transitions. 


“l Jenwary 1971 U. S. Merine Corps total was 21,3520. 


7 ——— 





TABLE II. MBASURES OF VARIATION 


Rank of VT NP £2/NP(I-P3) 
1 .014 3900 + 112 
Z none 4940 + 121 
3 TONS 3700 + 109 
4 SOL 7 - 2940 + 100 
5 vou 2660 + 96 
6 022 1860 + 82 


This completes the remarks on the steady state. In the 
following section, we shall consider the effect of an in- 
crease or decrease in the total number of officers in the 


SScemeaia Cvaluete transa@ent resuits 
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IV. ITIRANSTENT RESULTS 


Realistically, a graded personnel system which under- 
goes even occasional policy changes never actually reaches 
a steady state. Thus, analysis of the transient reaction 
We: the®system to a given policy Gifange is neceSisary if the 
eeerects of Such a change are to be well understood. | 

In (2) we defined eae as thegprobabality an ofbreen ts 
make at £ 1S j],.¢1ven he started inva at time zero. IThwvoushe 
out this paper we have assumed that Q(t) can be constructed 
from known data and can be considered as given. Even so, 
the direct determination of P(t) appears to be a difficult 
mask. Hence, it 15 desirable to develop a method to "“ac- 
Curateiy" approximate P(t). 

The relationship between P(t) and Q(t) was stated in 


(7) as 
P(t) = Hp)(t) + @ * P(t), 


and taking Laplace Stieltjes transforms and rewriting (11) 


we obtain 
P(s) = ) Q@(s)" p(s), 
n=0 : 
Ow tidt., arter anverting, 


prey = > 9) « H(t), (21) 
n=0 


where g(P) (4) is the nth fold convolution of oft) with it- 


self and ©) ey = I. 
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> 


: 


Smee =tnoewexpresslon 1o0r F(t) me (ai regmires tne 
Mateix convolution of @(t) with itself for all n = 1,2,3,... 
we seek simpler expressions which give bounds or approxima- 
mons tO P(t). 

tPEowwiersecertarniy an Upper bound for ce lise 
me SUDStitute P(t) = 7 in the repht hand side of (7) we can 


Write: 
Como nct) + G4 7 = Fea Hin), 


where H(t) is a matrix whose ith row contains elements H, (t). 
A second iteration, substituting the right hand side above 


moreet(t) yrelds 
Bet) <H,(t) tegen [2 yGe) mal coe 
< H(t) * [eS ey pe P(t) 1, 


ena byweoiuinuime this process, we have after ntl iterations 


nh 
p(t) < ) QO) « Hyct) + @™ 4 H(t). (22) 


As n>-o the term g™) nH (+) 0 SO that this @terative process 
G@@es converge to P(t) as defined in (21). 
Certainly, Bale 20, and an idemtieal iteratiame cal- 


CUlatroneshows that 


Th 
P(t) > » go ge aagek (23) 
k=0 


Thus, by combining Si and (23) we have 


ey) 


0 < P(t) ) glk) , H(t) < @ 1. (24) 


=0 





By choosing n large enough we can estimate P(t) by 


nN 
‘y gk) * H(t) 
k=0 


to any desired accuracy. 

In Part II a system of constant size N was considered 
and no changes in total size over time were investigated. 
bet us now comsider an increase in total size generated by 
increased inputs into state 1, i.e., the number of new re- 
cruits to the system exceeds the number departing from the 
system in a time period. This may be done on a one-time 
Deasis at a time we shall call zero, or by adding different 
amounts each year tor several successive years, starting at 
meee Zeno. “Filetmre 3 cdepicrs theme pe increase. 

Let n(t) = the number added at time t which are in 
excess Of those needed to replace departures from the system. 
Then N(t), the total system size at time t may be defined 


as 
t 


N(t) = N+ y n(u). (25) 


u=0 
This system is still "closed" and the concept of immediate 
replacement of departures applies for the entire N(t). If 
we assume the original system of N officers was in the steady 
State at time zero, the analysis of that group can be come 
ducted as discussed in Part III. To derive information 
concerning the entire system we can superimpose the tran- 
sient results of the added cohorts onto the o1 ginal system 


results. 


ao 





ADDED SIZE 


DUE TO 
INCREASE 
In(i) | | n(1) | 
[_n(0) | [nco) | 
—— —_ 
ORIGINAL 
SYSTEM 
Size 
N N N | N 
: as ic 
-1 0 1 2 


Hieine 5. oclze tinercace.. 


ree Y;(t) be the number of officers from those extra 
SOnorts added apwtime zerovand= after=tmateare in rank jaa 


jime t. Thus 


Y(t) =f), Yee 


Similarly define 5 05) and thie Vie CrOnmolen mero chesteotan 
system. Then 


ic 
EY; (t)] = )  n(t-v) Py; (u) (26) 


laid 
and thie weetor of expected numbers of e@fficers in the tomar 


system that are in the various ranks at time t is given by 
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Epewt)) ="EYX(=) )>) om (ey) 
NP eneceon 


NP, a ELY,(t) J 


(27) 


NP¢ + E[Y, (t) ] 


Hee SeCuCNat tomedwen late (27), De Clg is needed, and from 
(26) this requires knowledge of Pee u<t. These P44ft)'s 
can be approximated by (24). 

Mm discussed previously, the joint distribution oo: the 


Huper im cach rank at time t for a particular colort input 


1s multinomial. By assumption, each cohort is a MSMP which 
Pomindcpendciusor all other cohorts...) dhus 
; | 
VEY;(t)] = d a(t-u) P,;(u)(1-P;,(¥)), (28) 
and t 
cOV[Y, (t),¥,(t)] = - \ n(t-u) Py, Cu) Py, (u). (29) 
u=0 


Fommime the covariance matrix V[Y(t)] from (28) “and (23) 


me tollows that 


Ste = VEC Vee (30) 


To measure the effect of policy changes in the transient 
case, we must first translate these changes ino a new Q(t). 


Then, using (24}, we can estimate P(t). If policy changes 
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and an increase in overall size are to be effected simultan- 
eously, we can consider the entire system to be non-stationary 


and revise (26) to 


E{S(t)] = E[W(t)] + E[Y(t)] 


where 


6 
E[W,(t)] = % NP; Pgs i = 2y2 poate 
i=l 


Similarly, Ye replaces Pi in all calculations for V[X(t)] 
which is subsequently used in a revised (30). 

Tie Next agiemto be considéved@misea decrease in “oNvemael 
BPigeniilitlatea by awdecrease 1nweneanumber Of new officers 
Eieecting famk 1; that is, the number of officers leaving 
the system in a time period is greater than the number of 
replacements. This may be viewed as merely a diminished 
Pamlinile aput tosgcradually reduce size. .o7r as forced vac, 
Serc2on in the higher ranks with @ revised 90(t) that causes 
more "departures" to state 1, Whene™= Geplacement 15 m6 
Menger On an instantaneous one-to-one basis. <A major foxee 
reduction would undoubetdly be a combination of both policies. 

ket the size of the Officer Comps Prior to time zerogee 
N, and assume steady state conditions. Let N' < N be the 
desired new size. Then a system of size (N-N') must be 
emptied to reach this objective size. Figure 4 depicts thus 
peOCcess. 

At time zero we establish a subsystem of size N', whose 
members are chosen randomly such that the propertions of 


CeelWccws Mlecach rank are qual to tmose of tha tapoereen 


SZ 








aibglune sd, 


Ee 
>time 
c 


Size Reduction. 


Neoystem. These N' officers will proceed in the previously 


defined closed system model - in either the steady state or 


transient mode, as appropriate. 


However, the (N-N’) ether 


officers are placed in an open system, where state 0 is de- 


mime? ase OUlteOr the service." 


Distinguishing the 


Beeb y “, we Maye then @*=(t) for this open portion 


0 
Spl D 
ace) 
Cae) = Qzq(t) 
Q(t) 
ent) 
OPA (CE) 


FOr this Sap -system, HA Ct) 


| 


0 


Q12 
0 


0 


0 


3 F oT Ses ae 


and using the approximation (24) we can estima 


0 

(“e) 0 
Desi t) 

0 

0 

0 


0 


0 0 
0 0 
0 0 
Qz,(t) 0 
DS Sie 
0 0 
0 0 


t>o the limiting form of P*(t) is given by 


Dd 


. ps 


open sys- 


ao. 


(t) 0 


ect) 
0 


Thus, Ha (t) = ] 


(t).° As 





Paine) i 0 0 0 0 0 


t+ 





Tice tieme x mected TreMmainingenimecr In rank gq. 
BUR; tt), in the open system can De calculated using the 


Maisohenteparewstaces | through 6) of P*{7) 2 1..e., 


6 
= = t as 
ER; (t)] = ) (N-N') P, PE, (t). 
i=l 
AS Pe ad remk j E€mpties to zerGs ahus the expected 


number in the open subsystem, shown as the shaded area in 
Prepure 4, decreases to zero and the total objective size 


Nees peacned, | considering the system as a whole, we ave 
B(S;(t)] = g~ E{W,(t)] + E(R,(t)}; j = 1,2,...,6. 


TO conclude this section, we note that a logical mere 
Stcpelmmetie Utilization Of this model) as the construction 
Sec(tio tron actual policy and data. ~For cxample, consmden 
a simple development of OG eC) « Let 

O = aecrition rate from rank | due to death, disabilige 

diesmissal, Ctle., 

8 = opportunity for promotion. 

T= time in grade 1 until promotion atten < 


E@QSa:) 1. ~ 
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Then, 


Q17 0°) ss Qi y(t) fet recy = t<T 


1 (Teese a) (lee ie eee 


oP Gains t>EOS 
Q1 > (t) = QT, (t) =.0 t<T 
= B(1-a)! | tot 


Q(t) may then be substituted into (24) to obtain an 
estimate of P(t). Analysis of estimation errors involved 
and the number of convolutions necessary to achieve the 


desired accuracy can then be conducted. 
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